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Using a self-consistent Hartree-Fock approximation we investigate the relative stability of various stripe 
phases in the extended t-t'-U Hubbard model. One finds that a negative ratio of next- to nearest-neighbor 
hopping t' /t < expells holes from antiferromagnetic domains and reinforces the stripe order. Therefore 
the half-filled stripes not only accommodate holes but also redistribute them so that the kinetic energy is 
gained, and these stripes take over in the regime of t'/i — — 0.3 appropriate for YBa2Cu30B+i. 



Two-dimensional (2D) systems with strong electron correlations received a great deal of attention over 
the past twenty years. Such systems exhibit a subtle competition between rich variety of phases: an- 
tiferromagnetic (AF), ferromagnetic (FM), charge density wave (CDW), spin density wave (SDW), and 
superconducting phase. In particular, strong interplay of different order parameters in doped transition 
metal oxides leads to the experimentally observed phase separation involving holes ordered into domain 
walls (DWs) separated by AF domains 1 1 1. On the theoretical side, stripe phases have been predicted on 
the basis of Hartree-Fock (HF) calculations of doped AF clusters |°2l. One important feature of these results 
is that they indicate formation of the so-called^Weii stripes with density of one doped hole per one atom in 
a DW. However, neutron experiments on the Nd-codoped La2-a;Sr2;Cu04 (Nd-LSCO) cuprates reveal that 
the observed stripes are filled by one hole per two DWs, corresponding to half-filled stripes |3|. Unfor- 
tunately, the latter are only locally stable within the HF approximation 1 4 1 . Even though some additional 
Hamiltonian terms like the nearest-neighbor Coulomb repulsion V may slightly enhance their stability, 
filled stripes remain always better solutions f4l. Nonetheless, we will show below that under some circum- 
stances, half-filled vertical site-centered (HVSC) stripes may become energetically favored over the filled 
either diagonal site-centered (DSC) or vertical site-centered (VSC) ones even in the HF approximation. 

The starting point for the analysis of stripe structures is the effective 2D single-band Hubbard model. 



where the hopping ttj is t on the bonds connecting nearest-neighbor sites and t' for next-nearest-neighbor 
sites, while U is the on-site Coulomb interaction. The model is solved self-consistently within the HF 
method on a 16 x 16 cluster with periodic boundary conditions. For a representative doping level x = 1/8 
one obtains stable filled stripe structures with AF domains of width seven atoms. Their stability can be 
best understood from the band structure shown in Fig.^ In general, each filled DW induces the formation 
of two unoccupied bands, i.e., one for the f-spin and another one for the ^-spin, lying within the Mott- 
Hubbard gap, as illustrated in Fig.^a). Consequently, their special stability rests on a gap that opens in the 
symmetry broken state between the highest occupied state of the lower Hubbard band and the bottom of 
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Fig. 1 Band structure as a function of momentum ky, for a unit cell of width Nx = 16 in the Hubbard model with 
U/t = 5 for: (a) VSC stripe, doubled Brillouin zone; (b) HVSC stripe, quadrupled Brillouin zone. Black (gray) line 
corresponds to the bulk (mid-gap) bands, respectively, whereas the dashed line indicates the Fermi level. 



the mid-gap band. One may ask what happens now if one decreases the doping level so as to get half-filled 
stripes without a period quadrupling? Certainly, the twofold degenerate mid-gap states become quarter 
filled, as there is only half a hole per one DW. Hence, any stringent reason for such a symmetry breaking 
with the Fermi level in the middle of the lower mid-gap band is absent. Therefore, one needs to lower the 
symmetry by a period quadrupling so that the period on the stripe becomes four times the lattice constant 
l4l . This might be accomplished by introducing a SDW modulation along the DWs which leads to the 
band structure shown in Fig.^b), with a gap that opens up exactly at the Fermi energy. However, in spite 
of better optimizing the potential energy Eu, the HVSC stripe phase represents only a local minimum of 
energy, being less stable than both the filled VSC and DSC ones. Thus, guided by the observation that finite 
negative t' expels holes from the AF domains and reinforces the stripe order [5|, we investigate whether 
this mechanism suffices to stabilize the former structure in the t-t'-U model. 

In order to establish the role of t', we begin with energies of different phases att' = and t' = —0.3t 
(Table[0. The main effect is the energy gain of HVSC over DSC phase along (11) and (11) bonds with 
increasing \t'\. To understand it better we discuss Fig. |2l showing how the charge and spin configuration 
around the HVSC stripe is altered in the presence of increasing \t'\. In Fig. |2ja) we label all inequivalent 
bond kinetic energies along the (10) (A^;, . . . ,Da;), (01) (Ay, . . . ,Gy), (11) (A,. . . ,H), and (11) (in paren- 
thesis) directions. The most striking feature of the HVSC stripe phase is that the main kinetic energy gain 
is released not by transverse charge fluctuations, which is the case for the filled stripes, but by on-wall 
hopping processes on the bonds Cy (Fy) connecting sites with the FM (AF) coupling, respectively (see 
Table llJi. Another interesting property is that the system tries to regain some of the kinetic energy in the 

Table 1 Energies per site: ground-state energy -Etot, kinetic energy contributions for the bonds along (10) Ef, (01) 
E^, (11) SJT" and (11) directions, and potential energy Eu of the HVSC, VSC and DSC stripe phases in the 

t-t'-U model with U/t = 5 and s = 1 /8. VSC stripe is unstable at t' /t = -0.3. 





t'/t 


E-/t 


Ef/t 


E^-y/t 


E^'^+y/t 


Eu/t 


Etot/t 


HVSC 


0.0 


-0.5846 


-0.6585 


0.0000 


0.0000 


0.4627 


-0.7804 


VSC 


0.0 


-0.6753 


-0.6147 


0.0000 


0.0000 


0.4900 


-0.8000 


DSC 


0.0 


-0.6368 


-0.6368 


0.0000 


0.0000 


0.4696 


-0.8040 


DSC 


-0.3 


-0.6152 


-0.6152 


0.0000 


0.0257 


0.4379 


-0.7668 


HVSC 


-0.3 


-0.5786 


-0.6198 


-0.0080 


-0.0080 


0.4322 


-0.7678 
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Fig. 2 Rearrangement of the HVSC stripe order with increasing 
next-nearest-neighbor hopping: (a) t' /t = 0; (b) \t'\/t — 0.45. 
Circle sizes are proportional to doped hole density. Capital let- 
ters label bond kinetic energies along the (10) (Aj,, . . . ,Dx), (01) 
(Ay, . . . ,Gy), (11) (A,. . . ,H), and (11) (in parenthesis) directions, 
listed in Table|2| while + (— ) correspond to the bond kinetic energy 
gain (loss) with respect to the t' = case. Shadow circles indicate 
two crossing zigzag lines, see text (U/t=5, x — 1/8). 



Table 2 Bond kinetic energies along the 

(10) (A,,...,D,), (01) (Aj„...,G,), 

(11) and (11) (A,...,H) directions 
around the HVSC stripe shown in Fig.|2| 



t'/t 



bond 





-0.45 


A, 


-0.6081 


-0.6141 


B, 


-0.6652 


-0.5482 




-0.6277 


-0.6257 




-0.4377 


-0.4885 


Ay 


-0.6057 


-0.5822 


By 


-0.6003 


-0.5281 


c 


-0.8505 


-0.7234 


Dy 


-0.5759 


-0.5623 


F 

n^y 


-0.6515 


-0.5299 




-0.7206 


-0.6364 


Gy 


-0.6635 


-0.6046 


A 


0.0000 


-0.0053 


B 


0.0000 


0.0262 


C 


0.0000 


-0.1300 


D 


0.0000 


0.0128 


E 


0.0000 


0.0033 


F 


0.0000 


-0.0367 


G 


0.0000 


0.0164 


H 


0.0000 


0.0020 



transverse direction by developing a CDW on nearest-neighbor sites to the DW itself. Indeed, although 
the charge distribution along the stripe is uniform, the on-wall SDW causes the CDW on both sides of the 
DW. This clearly promotes the hopping between the AF sites (bonds B^) over the FM ones (bonds Dx). 
Furthermore, the CDW is strongly influenced by t', being first almost entirely quenched when t' ~ —0.15 
and then it is gradually restored [cf. Fig. EJb)]. 

These trends are fully consistent with the dynamical mean field theory studies of the effect of a single 
kink along the HVSC stripe (61 . It has been established that the energy cost of forming a kink increases 
slightly with increasing \t'\ up to \t'\/t < 0.1, whereas further increase leads to the opposite effect, and 
finally a wall with a kink becomes favored at \t'\/t ~ 0.3. Note, however, that although the overall shape 
of the brought back CDW is the same as of the initial one, the physical situation is fundamentally different 
- in the t' = Q limit, it is energetically advantageous for the system to equalize hole density between sites 
connected by the bonds (AF coupling), whereas in the large \t'\ limit, it tries instead to equalize hole 
density between sites connected by the bonds (FM coupling). Accordingly, such a charge redistribution 
results in two crossing zigzag lines, which facilitate hole propagation. Indeed, each path consists of a -C-F- 
bond pattern along the (11) direction and analogous sequence in the (ll) direction. Even though a moving 
hole gains kinetic energy on both type of bonds, a more significant gain is achieved on the FM bond C, 
being the driving force in the formation of zigzag pattern. It is also supported by a systematically growing 
kinetic energy gain with increasing \t' \ on the FM bond connecting the zigzag paths (cf. Table|2j. 

The effect of increasing \t'\ on the kinetic energy per diagonal bond E'f of the HVSC stripe phase 
illustrates Fig.|3ja). One observes that in the small \t'\ regime next-nearest-neighbor hopping processes 
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Fig. 3 Average kinetic energy Efi per diagonal bond as a function of increasing next-nearest-neighbor hopping \t'\ 
for: (a) half-filled vertical site-centered (HVSC) stripe phase; (b) filled diagonal site-centered (DSC) stripe phase, as 
obtained for U/t — 4, 5, 6, 7 (solid, dashed, long-dashed, dot-dashed line). Panel (c) shows the phase boundary 
between the DSC stripe phase and the HVSC one as obtained for doping x — 1/8. 



cost energy and therefore they are suppressed by quenching the CDW order on sites next to the DW. 
However, towards larger \t'\, the kinetic energy associated with next-nearest-neighbor hopping becomes 
negative despite the negative sign of t'. Consequently, the system develops a new CDW so as to optimize 
the Ef' gain. In contrast, Et' of the filled DSC stripe phase remains positive even in the large \t'\ = OAt 
regime, as shown in Fig.|3jb). Moreover, increasing Coulomb repulsion U yields in this case only a small 
suppression of Ef, whereas the same increase of U is reflected in much larger Et' gain of the HVSC stripe 
phases. This, together with a better on-site energy Ejj, explains a broader region of stability of the HVSC 
stripe phases in the strongly coiTelated regime, as illustrated in Fig. |3jc). We interprete their enhanced 
stability in the small U regime as following from the melting of both stripe phases. Indeed, increased 
mobility of the holes released from the DSC stripes results in a faster, as compared to the HVSC case, 
enhancement of the Ef energy cost upon increasing \t'\. 

In conclusion, we have found that the next-nearest-neighbor hopping t' plays an important role in af- 
fecting the relative stability between filled and half-filled stripe phases. However, the established value 
\t'\/t ~ 0.3 opening a window for the stability of the HVSC stripe phase, is excessively larger than 
t' /t —0.1 of LSCO compounds Q. Interestingly, t' /t ~ —0.3 corresponds to the value appropriate 
for YBa2Cu306+5 (YBCO) and indeed a half-filled vertical stripe phase has been reported in this com- 
pound for S = 0.35 1 8 1. Therefore, our result calls for a further experimental characterization of stripes in 
YBCO. Finally, it may be expected that a better treatment of electron coiTelations would stabilize further 
the half-filled stripe phases (9) for realistic values of t'. 
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